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Recently, there have been suggestions that the apparent accelerated expansion of the universe is 
not caused by repulsive gravitation due to dark energy, but is rather a result of inhomogeneities 
in the distribution of matter. In this work, we investigate the behaviour of a dust dominated 
inhomogeneous Lemaitre-Tolman-Bondi universe model, and confront it with various astrophysical 
observations. We find that such a model can easily explain the observed luminosity distance-redshift 
relation of supernovae without the need for dark energy, when the inhomogeneity is in the form of 
an underdense bubble centered near the observer. With the additional assumption that the universe 
outside the bubble is approximately described by a homogeneous Einstein-de Sitter model, we find 
that the position of the first CMB peak can be made to match the WMAP observations. Whether 
or not it is possible to reproduce the entire CMB angular power spectrum in an inhomogeneous 
model without dark energy, is still an open question. 
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I. INTRODUCTION 


The first indications that the universe is presently in a 
state of accelerated expansion were given by J.-E. Sol- 
heim as far back as in 1966 Q. Using the observed 
luminosity of several cluster galaxies he found that the 
model giving the best fit to the data was one with a 
non-vanishing cosmological constant and negative decel¬ 
eration parameter. It is, however, only after the more 
recent observations of the luminosity of supernovae of 
type la (SNIa) that this claim has grown in popularity. 
The first SNIa observations supporting this claim were 
those of Riess et al. in 1998 Q and Perlmutter et al. in 
1999 Q- Since then, more recent observations of smer- 
novae seem to strengthen this claim even further . 

Other independent observations that appear to favour 
the picture of a universe in a phase of accelerating expan¬ 
sion, are the measurements of the anisotropies in cosmic 
microwave bacteround (CMB) temperature 0 and the 
galaxy surveys [g • With these observations in mind, the 
current period of accelerated expansion seems to be well- 
established. The physical mechanism that drives this ac¬ 
celerated expansion is, however, still an open question. It 
is usually ascribed to an exotic energy component dubbed 
dark energy, whose nature remains a mystery. 

Recently, there have been several papers discussing the 
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possibility that the apparent accelerated expansion of the 
universe is not caused by this mysterious dark energy, 
but rather by inhomogeneities in the distribution of mat¬ 
ter. Most of these papers look at the backreaction effects 
arising from perturbing homogeneous models, and try to 
explain the accelerated expansion as corrections to the 
zeroth order evolution from the higher-order, inhomoge¬ 
neous terms (see e.g. 0,0,O,0,|3Q,0)- However, 
several pap ers criticizing some of this work have appeared 
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Another approach is to look at inhomogeneities of a 
larger scale in the form of underdense bubbles. The ba¬ 
sic idea behind this line of explanation is that we live in 
an underdense region of the universe, and the evolution 
of this underdensity is what we perceive as an acceler¬ 
ated expansion. An analysis of early supernova data by 
Zehavi et al. gave the first indications that there might 
indeed exist such an underdense bubble centered near us 



Specific models that give rise to such underdensities have 
been studied prev iousl y in the form of a local homo¬ 
geneous void |22l 111 In these works both the 

underdensity and the region outside it are assumed to 
be perfectly homogeneous Friedmann-Robertson-Walker 
(FRW) models with a singular mass shell separating the 
two regions. The inhomogeneity manifests itself as a dis¬ 
continuous jump at the location of the mass shell. 

In this article, we wish to investigate more realistic mod¬ 
els where there is a continuous transition between the 
inner underdensity and the outer regions. Therefore we 
consider an isotropic but inhomogeneous dust dominated 
universe model, where the inhomogeneity is spherically 
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symmetric. The model can then be described within the 
Lemaitre-Tolman-Bondi (LTB l cl ass of spherically sym¬ 
metric universe models j2^|^|^ . To make contact with 
the ordinary FRW models, we assume that the universe 
is homogeneous except for an isotropic inhomogeneity of 
limited spatial extension, where the transition between 
these two regions is continuous. 

In a homogeneous universe, it is possible to infer the 
time evolution of the cosmic expansion from observa¬ 
tions along the past light cone, since the expansion rate 
is a function of time only. In the inhomogeneous case, 
however, the expansion rate varies both with time and 
space. Therefore, if the expansion rates inferred from 
observations of supernovae are larger for low redshifts 
than higher redshifts, this must be attributed to cosmic 
acceleration in a homogeneous universe, whereas in our 
case it can simply be the consequence of a spatial vari¬ 
ation, with the expansion rate being larger closer to us. 
As shown in [^, this results in an expression for the 
luminosity distance-redshift relation where the inhomo¬ 
geneity mimics the role of the cosmological constant in 
homogeneous models. 

However, the supernova observations are not the only 
data that support the claim of an accelerating expansion. 
As mentioned above, CMB observations also seem to lend 
support to this claim. Therefore, in order for our model 
to be considered realistic, it should also be able to explain 
the observed CMB temperature power spectrum. We will 
not attempt to reproduce the whole CMB temperature 
spectrum for our inhomogeneous model in this paper. 
For simplicity, we will limit ourselves to the location of 
the first acoustic peak. As we will show in section Hill 
it is possible to obtain a very good match to both the 
supernova data and the location of the first acoustic peak 
simultaneously. In fact, the match to the supernova data 
is better than for the ACDM model. 

The observed isotropy of the CMB radiation implies that 
we must be located close to the center of the inhomo¬ 
geneity. According to this picture, we are positioned at a 
rather special place in the universe. On the other hand, 
this model has the attractive feature that there is no need 
for dark energy. Also the model is sufficiently simple so 
that it can be solved exactly. It is therefore a good toy 
model for testing the ideas of inhomogeneities as a solu¬ 
tion to the mystery of the dark energy. 

The structure of this paper is as follows: First, we will 
present our model in section im parameterized by two 
functions a(r) and /3(r) related to the distribution of 
matter and spatial curvature, respectively. Still in sec¬ 
tion QH we present the formalism needed in order to ob¬ 
tain the luminosity distance-redshift relation for spheri¬ 
cally symmetric, inhomogeneous models. In section uni 
we discuss the physics behind the first peak of the CMB 
spectrum and define a shift parameter that quantifies the 
deviation of the location of this peak relative to that of 


the concordance ACDM model. In section Hvl we present 
the results from the confrontation of our model with the 
physical tests presented in the preceding sections, and 
discuss briefly the possibility of using the recently de¬ 
tected baryon oscillations in the matter power spectrum 
to constrain the model even further. Finally, in sectionIVI 
we summarize our work. 


II. SPHERICALLY SYMMETRIC, 
INHOMOGENEOUS UNIVERSE MODELS 


The line element for a general, spherically symmetric, 
inhomogeneous universe model may be written 

ds^ = —dt^ + A^(r, t)dr'^ + R^{r, t)d^^ . (1) 


The Einstein equations are 

G = R^u ~ ( 2 ) 

where k = SttG and the energy-momentum tensor is as¬ 
sumed to be = diag(p, 0, 0, 0), i.e. containing matter 
only. 


Solving the equation Gqi = 0 gives 


X{r,t) 


R'{r,t) 


( 3 ) 


where (3(r) is an arbitrary function of r. Throughout this 
paper, we will use a ' = d/dr to denote differentiation 
with respect to r and ' = d/dt for differentiation with 
respect to t. 

The Einstein equations for the dust dominated Lemaitre- 
Tolman-Bondi universe models take the form 

(4) 

-67Li(7^ + 27Li-2^-27L||Hj. + ^ = -«p(5) 

where = R/R, iLy = R'/R' and q±_ = —RR/R^. 
Adding Eqs. 0 and ®, we obtain 

2RR + R^ =(d (6) 


Integration of this equation leads to 


Hi 


(3 a 


( 7 ) 


where a is a function of r. 


Hence, the dynamical effects of (3 and a are similar to 
those of curvature and dust, respectively. 
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Differentiating Eq. 0 with respect to r and inserting the 
result into Eq. we obtain the density distribution as 

= isif <«> 


Substituting Eqs. 0 and o into the expression for the 
deceleration parameter yields 

1 a 1 a , . 

Obviously, this quantity is non-negative (since a > 0) 
and equal to the usual value q± = 0.5 for a spatially flat, 
dust dominated universe. Thus, an inhomogeneous, dust 
dominated universe cannot be accelerating in the sense 
of having a negative deceleration parameter. 


Since we are interested in the late time behaviour of this 
model, we will define f = 0 as the time when photons 
decoupled from matter, i.e. the time of last scattering. 
Furthermore, we define R{r,t = 0) = Ro{r) and intro¬ 
duce a conformal time ry by = Rdrj. Integrating 

Eqs. 0 and 0 with /3 > 0 yields 


R 




a 

w 


(coshry — 1) 


-|-i?o 


cosh 7y 


’ a + /3i?o . , 


— (sinh7y-ry) 


-|-i?0 


sinh 7y - 


/ o -t- j3R(j 

pRo 


(coshTy — 1) 


( 10 ) 


( 11 ) 


which is an exact solution of Einstein’s equations for this 
class of models. 


The ordinary dust dominated solution for a universe with 
negative spatial curvature is found by choosing 

a = Hgn^or^, P = Hgidkor'^ and i?o = 0 (12) 

yielding 

R = ^ 7 f^(coshry-1) = ra(ry) (13) 

^ (14) 

2Ho 

and Kp = . Here, a(?y) is recognized as the 

scale factor in the FRW model, while Hmo and 0.^0 are 
the matter and curvature density today, respectively. 


Since we are interested in studying a universe with an un¬ 
derdensity at the center, we choose the a and P functions 
so that they interpolate between two such homogeneous 
solutions: 


a{r) 

P{r) 


— 


Qfo - Aa ( - 


= /3o - A/3 ( ^ 


- tanh 


— tanh 
2 


2Ar ) 
r-ro \ 
2Ar J 


(15) 

(16) 


Here, JI±^o is the value of the transverse Hubble parame¬ 
ter in the outer homogeneous region today, while ap and 
Pp are given by the matter and curvature density in this 
region, respectively. Furthermore, Aa and A/3 specify 
the differences in the parameters between the two re¬ 
gions, and rp and Ar specify the position and width of 
the transition. 

The function Rp{r) can be chosen freely by a suitable 
choice of coordinates r (if the universe has a finite size 
at f = 0). To match our solution to a homogeneous 
FRW solution in the outer region, we choose Rp = a^:r, 
where a* is the scale factor of the homogeneous model at 
recombination. 

To relate the a and /3 functions to observable quantities, 
we define relative matter and curvature densities from 
the generalized Friedmann equation 0 as 




m 




Kp 

Hl + 2H^ 

1 Hm 


(17) 

(18) 


Note that for the homogeneous case, with H± = 
these expressions coincide with the usual definitions. 

Furthermore, we need to find the luminosity distance- 
redshift relation in this model for comparison with su¬ 
pernova observations. The photons arriving at r = 0 
today (defined as t = tp) follow a path t = t{r) given by 


di R'{r,t) 
dr + P 


(19) 


with f(0) = to- Following Iguchi et al. [^, we find the 
redshift z = z{r) of these photons from 


dz 

dr 


= {z + l) 


■\/l + /I 


( 20 ) 


with the initial condition z(0) = 0. The position of the 
last scattering surface (i.e. the position of the CMB pho¬ 
tons that we observe today, at the time of last scat¬ 
tering) is given by t(r*) = 0, and we define tp by 
z{rif) = Zt. ~ 1100. An accurate formula for z* in terms 
of the matter contents of the universe has been given by 
Hu and Sugiyama 

The luminosity distance is then given by the usual ex¬ 
pression 


dkiz) = (1 -H zfR{r, i) (21) 

and the angular diameter distance is 
dA{z) = Rir,i) 


( 22 ) 
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III. THE COSMIC MICROWAVE 
BACKGROUND 


To confront our model with observations of the CMB, 
we would in principle need to study perturbations in an 
inhomogeneous universe. However, since our model is 
homogeneous outside a limited region at the center, we 
will assume that the evolution of perturbations is iden¬ 
tical to that in a homogeneous universe until the time 
of last scattering. This means that we can use the ordi¬ 
nary results for the scale of the acoustic oscillations at 
the last scattering surface. On the other hand, the an¬ 
gular diameter distance, which converts this scale to a 
corresponding angle on the sky, is sensitive to the inho¬ 
mogeneity at the center. As a simple test we will use 
the position of the first peak in the CMB spectrum to 
constrain our inhomogeneous models. 

The position of the m-th Doppler peak in the CMB spec¬ 
trum can be written as 

^m)^A 

where I a is the acoustic scale and Lpm is a small shift 
mainly due to the projection of the three-dimensional 
temperature power spectrum onto a two-dimensional an¬ 
gular power spectrum. 

The acoustic scale is given by 


where (Ia is the angular diameter distance to the last scat¬ 
tering surface and is the sound horizon at recombina¬ 
tion. In a standard FRW cosmology, these two quantities 
are approximately given by (in comoving coordinates) 


where oi and 02 are given by 

oi = 0.286-b 0.626^6 (29) 

02 = 0.1786- 6.308a;b-bl74.9a;^-1168a;^ (30) 

To reduce the effect of the approximations made in the 
above formulae, we will introduce a shift parameter that 
measures the position of the first Doppler peak for a given 
model relative to the concordance ACDM model. That 
is, we define 

5 = (31) 

where is the peak position for the current con¬ 

cordance model, with uJk = 0, uJm = 0.135, cut = 0.0224, 
oj^ = 4.2 • 10“® and lua — 0.368. To be in agreement with 
the WMAP observations, the shift parameter should be 
within the range S = 1.00 ± 0.01. In fact, the relative 
error in the peak position from the WMAP data is 
0.8/220.1 ~ 4-10“^. However, the approximations made 
in the formula in Eq. m are probably of the same or¬ 
der of magnitude. Therefore it is safe to say that models 
with |5 — 1| > 0.01 are ruled out, whereas models with S 
within a percent of the ACDM value are probably worth 
a closer look. After all, there is still a long way to go 
from the correct position of one peak to a perfect match 
with the entire CMB angular power spectrum. 

In addition to the correct peak position, cvh should be 
within the range predicted by Big Bang nucleosynthesis 
1^. ujb = 0.020 ± 0.002. For simplicity, we will use the 
best-fit value given by the WMAP team 0, Wb = 0.0224. 

Inserting the value for , Eq. (mil becomes 

5 = 0.01419(1-^51)— (32) 


r - r 
" /o ait') 


and 


dA = 




-rSk 


feOl 




dz 


feol 


H{z) 


(25) 


(26) 


where Cg is the sound speed of the baryon-photon plasma 
prior to recombination and z* = 1 /a, — 1 is the redshift 
at the time of recombination. The function Sk depends 
on the spatial curvature and is defined as 


Skix) 


sin a; , flk < 0 

X , Dfc — 0 

sinha; , 0 ^ > 0 


(27) 


A fitting formula for the dependence of ipi on Ub and Urm 
where uji = Q,ih? is the density of energy component i, 
can be found in Ref. . The formula can be written as 


ipi = ai 


UJ 


7 


02 




(28) 


As an example, an Einstein-de Sitter (EdS) model with 
Um = 1, Hubble parameter h — 0.71 and ojb = 0.0224 has 
S = 0.916, whereas the same model with h = 0.51 has 
S = 0.998. In comparison, CMBFAST yields the val¬ 
ues S — 0.914 and S = 0.998 for these two models. As we 
can see, the formalism accurately describes the position 
of the first peak. The CMB temperature power spectra of 
these three models are plotted in figure ^ Note that the 
second and third peaks will also have approximately cor¬ 
rect positions when the shift parameter is close to 1 , since 
the in Eq. (HI are relatively small. For instance, 
the EdS model with h = 0.51 has I 2 = 530 and Z 3 = 784, 
whereas the best-fit values from are I 2 = 529 and 
?3 = 781. Also note that the relative error in the location 
of the second and third peak are larger than for the first 
peak, at around 3%. 

The relative heights of the peaks, on the other hand, 
have a more complicated dependence on the parameters 
of the model, see e.g. |^. We will postpone discussing 
these features of the CMB spectrum until we have a bet¬ 
ter understanding of the evolution of perturbations in an 
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FIG. 1: CMB angular power spectra for the ACDM model 
and two Einstein-de Sitter models, normalized to the height 
of the first Doppler peak 

inhomogeneous model. Note, however, that it is possible 
to make matter-dominated homogeneous models that fit 
the observed CMB spectrum, see e.g. |^ . 

In our case, we must use expression for the angular 
diameter distance to the last scattering surface when we 
calculate the shift parameter in Eq. E3)- On the other 
hand, we can still use the expression for the sound horizon 
as defined in the homogeneous case in Eq. , since our 
model is assumed to be homogeneous close to the last 
scattering surface. 


IV. RESULTS 


When going from a homogeneous to an inhomogeneous 
universe model, the parameters describing the model (tUm 
and LOk) become functions of r. This means that we in¬ 
troduce, in principle, an infinite number of new degrees 
of freedom. However, for the purpose of studying the 
possibility of explaining the current observations with¬ 
out introducing dark energy into the model, we have re¬ 
stricted ourselves to a very simple “toy model”: An un- 
derdense region close to us, surrounded by a flat, matter 
dominated universe. This means that we must choose 
Qfo = 1 and /3o = 0. Furthermore, we put Aa = —A/3. 
This leaves four parameters, Aa, ro, Ar and the physical 
Hubble parameter at the origin, H||(0,to) = lOOfe 
to be fitted to the observations. 

Let us first focus on the two main observations: The su¬ 
pernova Hubble diagram and CMB angular power spec¬ 
trum. A good fit to the supernova data requires the Hub¬ 
ble parameter inside the underdensity, /lin to be around 



FIG. 2: The spatial variation of the Hubble-parameters at 
t = to 


/ijn ~ 0.65. On the other hand, a good fit to the CMB 
spectrum for a flat matter dominated model requires 
the Hubble parameter outside the underdensity to be 
/lout — 0.5. This more or less determines the two pa¬ 
rameters Act and h. 

Next, the shape of the transition between the underdense 
and the homogeneous region is specified by rg and Ar. 
These values dictate the redshift-magnitude relationship, 
and must be chosen to fit the supernova Hubble diagram. 
There are lots of choices for the parameters that give 
a very good fit to both the supernovae and the posi¬ 
tion of the first acoustic peak in the temperature power 
spectrum. However, we want the underdensity in our 
model to be such that the matter density is compati¬ 
ble with the current model independent observations of 
rirno- An excellent candidate for such observations is the 
mass-to-light ratio measurements made by the 2dF team 
[^. These yield ^mo = 0.24 ± 0.05 from observations 
of galaxies with redshifts z < 0.12. We will therefore 
choose the free parameters such that the mass density 
parameter at the origin is within this range in addition 
to giving a good fit to the supernova measurements and 
the CMB peak. The model which we adopt as our “stan¬ 
dard model” gives a matter density at the center of the 
underdensity of flmo = 0.20. A plot of the spatial vari¬ 
ation today of the Hubble parameters of our standard 
model is given in figure|21 Furthermore, a plot of the dis¬ 
tance modulus of this model together with the supernova 
observations can found in figure 13 . 

Note that the y^-value for our model is = 176.5, when 
compared to the “gold” dataset of Riess et al. . This 
is sligh^ better than that of the concordance ACDM 
model 0, x\cDM = 178. 

The spatial shapes of the underdensity at the initial time 
and today are plotted in figure^las functions of the phys- 
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FIG. 3: Distance modulus vs. redshift for our standard model 
together with supernova observations 



FIG. 4: The evolution of the underdensity in our standard 
model 


ical distance today. This illustrates the time evolution of 
the underdensity. As we can see, the shape stays almost 
constant. This is due to the Hubble parameters H|| and 
being roughly constant in space. 

Although the matter distribution is clearly inhomoge¬ 
neous close to the observer, we wish to point out that 
this does not necessarily contradict the data from the 
galaxy surveys. One often hears the claim that these 
surveys show the local distribution of matter to be ho¬ 
mogeneous. However, it is probably more correct to say 
that they are shown to be compatible with a homoge¬ 
neous universe rather than actually proving it. The key 
point here is that in order to determine for example the 
number counts of galaxy clusters, one needs to make an 


assumption about galaxy evolution and how likely it is 
to observe a galaxy with a certain luminosity at a cer¬ 
tain redshift. As pointed out in e.g. one usually 

assumes a homogeneous universe in order to deduce the 
effects of source evolution. Therefore, using this deduced 
evolution to claim observed homogeneity in the number 
counts amounts only to circular arg umentation. Further¬ 
more, it is explicitly shown in |39l| that given any LTB 
model it is always possible to find a source evolution that 
agrees with the observed number counts. 

The inhomogeneity at the center gives only a minor 
change of the angular diameter distance to the last scat¬ 
tering surface. In fact, our model has (Ia = 10.4 Mpc, 
which is the same value as we find for the Einstein-de 
Sitter model with h = 0.51. (Note that these values 
are physical, not comoving distances). Our model, with 
Hubble parameter h ~ 0.51 in the homogeneous region, 
thus yields a CMB angular power spectrum very similar 
to the one plotted in figure Q] at least for large I values. 
Using the formula E3), we find S = 1.006, i.e. an almost 
perfect match for the position of the first Doppler peak. 
For smaller I values, the CMB pattern will be affected 
by our position relative to the center of the underden¬ 
sity. This has been studied in the previously mentioned 
void model of Tomita , who concluded that relatively 
large displacements from the center of the underdensity 
were fully consistent with the observed CMB dipole and 
quadrupole. Furthermore, J. Moffat ^ argues that such 
a displacement could even explain the detected alignment 
of the CMB quadru- and octopole ^3 ■ 

A rough estimate of the apparent peculiar velocity for an 
off-center observer is 

Vp ^ (h||_in - /i||,out) lo ■ lOOkm/s (33) 

where Iq is the distance from the observer to the center, 
measured in Mpc. If we for instance require that Vp must 
be less than the estimated peculiar velocity of the local 
group which is of the order of 600km/s, this means 
that the observer must be within 40 Mpc from the center 
of the inhomogeneity. Even stronger constraints might be 
obtained by considering the peculiar velocities of nearby 
clusters, see e.g. [d^. 

Recently, Eisenstein et al. announced the detection of 
baryon oscillations in the SDSS galaxy power spectrum 
[ 43 . This represents additional, independent data that 
can be used to constrain our model even further. The 
physical length scale associated with these oscillations is 
set by the sound horizon at recombination. Measuring 
how large this length scale appears at some redshift in 
the galaxy power spectrum allows us to constrain the 
time evolution of the universe from recombination to the 
time corresponding to this redshift |^l4l . l4^l4^l4^l43| . 

A length scale quoted by Eisenstein et al. is the ratio of 
the effective distance to the chosen redshift in the galaxy 
survey to the angular diameter distance to the last scat- 
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Description 

Symbol 

Value 

Density contrast parameter 

Ao 

0.90 

Transition point 


1.35 Gpc 

Transition width 

Ar/ro 

0.40 

Fit to supernovae 

Xsw 

176.5 

Position of first CMB peak 

S 

1.006 

Age of the universe 

to 

12.8 Gyr 

Relative density inside underdensity 

^m,in 

0.20 

Relative density outside underdensity 

,out 

1.00 

Hubble parameter inside underdensity 

hin 

0.65 

Hubble parameter outside underdensity 

hout 

0.51 

Physical distance to last scattering surface 

Dlss 

11.3 Gpc 

Length scale of baryon oscillation from SDSS 

Ro.35 

107.1 


ated expansion of the universe without the need to in¬ 
troduce dark energy. Inspired by the recent discussions 
about the possibility of explaining the apparent acceler¬ 
ation by inhomogeneities in the matter distribution, we 
have studied a model where the observer is assumed to 
be situated near the center of an underdense bubble in a 
flat, matter dominated universe. If this model is realistic, 
we live in a perturbed Einstein-de Sitter universe within 
130 million light years from the center of an underdensity 
that extends about 5 billion light years outwards. Under 
the assumption of spherical symmetry this universe is 
described by the Lemaitre-Tolman-Bondi space-time. 


TABLE I: The parameters and features of our adopted stan¬ 
dard inhomogeneous model 


tering surface, 


^ 0.35 — 


dvi^sdss^ 

dA{z^) 


(34) 


where Zsdss = 0.35. The effective distance dy is defined 
in Eq. (2) in Ref. as a mix of radial and angular 
distance, to take into account that these scale differently. 
The value they measure for this ratio is i?o .35 = 79.0 ±2.9. 
Note that this value differs from that quoted in The 
reason for this is that we’ve chosen to give the distances 
dv and d^ in physical coordinates, while Eisenstein et al. 
quote them as comoving. 


Calculating this ratio for our model we find the value 
= 107.1. Comparing this value with that quoted 
by Eisenstein et al., one might be tempted to claim that 
the model is ruled out. However, in order to say some¬ 
thing conclusive using this constraint, we need to be sure 
that the “measured” value of i?o .35 is model indepen¬ 
dent. But when the authors derived this constraint they 
assumed a ACDM model. This makes it a little unclear 
how to use this constraint for non-ACDM models, or, in¬ 
deed, whether it is even possible to use it for such models. 
Ideally, one would need to repeat the analysis of Eisen¬ 
stein et al. assuming our inhomogeneous model as base 
model. We will therefore be careful not to rule out the 
model based on this parameter alone. 


The main features of our standard model are summarized 
in tabled Note that the age of the universe is 12.8Gyr 
in our model. This is significantly less than the value for 
the concordance ACDM model, 13.7 Gyr, but it is still 
in agreement with observations of globular clusters [i^ . 
which put a lower limit of 11.2 Gyr on the age of the 
universe. 


V. SUMMARY AND CONCLUSIONS 


The main goal of this paper has been to present a simple 
model with the ability to explain the apparent acceler- 


The two main observations we sought to explain were the 
luminosity distance-redshift relation inferred from SNIa 
observations and the GMB temperature power spectrum. 
These two sets of observations are made at opposite ends 
of the redshift spectrum, respectively low redshifts for 
the supernovae and high redshifts for the GMB. The fact 
that our model is inhomogeneous allows us therefore to 
choose the geometry and matter distribution such that 
the physical conditions are favourable for explaining the 
SNIa at low redshift while they at the same time are 
favourable for explaining the GMB at high redshifts. We 
find that a very good fit to the supernova data is obtained 
if we allow the transverse Hubble parameter to decrease 
with the distance from the observer. On the other hand, 
we get a good fit to the location of the first peak of the 
CMB power spectrum if we assume the universe to be flat 
with a value of 0.51 for the Hubble parameter outside the 
inhomogeneity. Interpolating between these two limiting 
behaviours we get a good fit to both the supernova data 
and the location of the first peak. 

Our model yields a better fit to the Riess data set of su¬ 
pernovae than the concordance AGDM model. However, 
for the GMB fit we tested only for the location of the first 
peak. Although the model yields a good fit to this, it does 
not necessarily mean that it matches the whole GMB 
spectrum. Indeed, since the physics responsible for the 
acoustic peaks is determined by the pre-recombination 
era, we would expect the peaks to look more or less the 
same as for a flat, homogeneous model with h = 0.51. 
This suggests that our model might fail to explain the 
third peak. Furthermore, the model does not appear to 
be able to explain the observed length scale of the baryon 
oscillations in the SDSS matter power spectrum either, 
although one may question whether the data quoted by 
the SDSS team can be used directly to test our model. 

The most powerful way to rule out inhomogeneous uni¬ 
verse models would be to do a full analysis of the evolu¬ 
tion of perturbations in these models. In that way, one 
could confront the model with both the full GMB angular 
power spectrum and the matter power spectrum. Only 
after such an analysis is carried out can one say whether 
our model is ruled out or if it is a viable alternative to 
dark energy. 
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